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The elastic and inelastic scattering of protons from mirror nuclei ^**Ne and ^**0 are studied 
in a folding model approach. For comparison, two different effective interactions are folded with 
Hartree-Fock densities to obtain the nuclear interaction potentials. Both of them provide equivalent 
descriptions to the data and the deformation parameters extracted from inelastic scattering are 
reasonable. The density dependence parameters obtained from nuclear matter calculations and 
used for present analysis also provide a good estimate for the nuclear mean free path. The present 
formalism unifies radioactivity, nuclear matter and nuclear scattering. 
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O '■ I. INTRODUCTION 

D ■ 

' Proton scattering has been widely used as a means to study both collective and microscopic aspects of nuclear 
Q\ ' structure [1,2]. The study is consistent only if a well-defined effective nucleon-nucleon (NN) interaction is applied 
in the analysis. Also, with the advent of radioactive nuclear beams there is constant enhancement of our knowledge 
frontiers on the structure and reaction dynamics of the known stable nuclei as well as their less known unstable 
counterparts [3,4]. This rapidly developing field provides a testing ground for different nuclear reaction theories and 
effective interactions. Scattering involving ^*Ne and are interesting because not only are they mirror nuclei, but 
also ^^O is a stable nucleus while ^*Ne is its radioactive counterpart. 

In this work, proton scattering on ^^Ne and has been studied at low energies (<100 MeV/A) [5,6] in a folding 
\ model approach. The folding model is well known as a powerful tool for analyzing nucleus- nucleus scattering data 
• at relatively low incident energies [4,7-9]. It directly links the density profile of the nucleus with the scattering 
[ cross sections and is thus very appropriate for studying nuclei, especially those with exotic matter distributions. A 
semi-microscopic analysis in the optical model (OM) framework is carried out. In the DWBA calculations of nuclear 
excitation, with transferred angular momentum I, the form factors are obtained by taking the derivatives of the 
potentials used. 
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II. THEORETICAL FORMULATION 

The nucleon-nucleus potential can be obtained by single folding the density distribution of the nucleus with the 
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^ , nucleon-nucleon effective interaction [10] as. 



U{ri)^ J P2{ri)vNN{\rl -r-i\)d^r-i (1) 

where P2{f^2) is density of the nucleus at f2 and unn is the effective interaction between two nucleons at the sites f\ and 
r2 . Two different forms of effective interactions have been employed in this work. We perform a comparative study 
between the Modified Seyler-Blanchard (SBM) and density dependent M3Y (DDM3Y) effective NN interactions. 

The finite range, density, momentum and isospin dependent effective interaction SBM has different strengths for 
pp (or nn) and pn interactions and its form is [11], 

v{r = |rl - rl|,p, p) = -Ci.J—^\\ -^-d'{p,+ p,)-] (2) 

r/a 0^ 
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where, the subscripts T and 'u' refer to hke-pair (nn or pp) and unUke-pair (np) interactions, respectively. Here 'a' 
is the range of the two-body interaction, '6' is a measure of the strength of repulsion with relative momentiun 'p\ 
while 'd' and 'n' are two parameters determining the strength of density dependence. pi{fi) and P2ii^) arc densities 
at the sites of the two interacting nucleons. The values of parameters n, Ci, C^, a, b, d are given in Table 1. These 
constants are found to reproduce the bulk properties of nuclear matter and of finite nuclei [11,12] and are known also 
to explain the p + 4,6,8jjj, 6,7,9,11^^ scattering data successfully [2,4,13-16]. The parameters are determined without 
exchange effects and thus they contain the effect indirectly though in a very approximate way. 
The finite range M3Y effective interaction v{r) appearing in Eqn. 1 is given by [17] 

vir) = 7999 2134 (3) 

^ ' 4r 2.5r ^ ' 

This interaction is based upon a realistic G-matrix, which was constructed in an oscillator representation. Effectively 
it is an average over a range of nuclear densities as well as energies and therefore the M3Y has no explicit density or 
energy dependence. The only energy dependent effect that arises from its use is a rather weak one contained in an 
approximate treatment of single-miclcon knock-on exchange. The density and energy averages are adequate for the 
real part of the optical potential for heavy ions at lower energies. Although, it is important to consider the density 
and energy dependence explicitly for scattering at higher energies, where the effects of a nuclear rainbow are seen and 
hence the scattering becomes sensitive; to the potential at small radii. Such cases were studied introducing suitable 
and semirealistic explicit density dependence [18,19] into the M3Y interaction which was then called the DDM3Y and 
was very successful for interpreting consistently the high energy elastic a and heavy ion scattering data. The present 
calculations use the density dependent M3Y effective interaction (DDM3Y) supplemented by a zero-range pseudo 
potential. In DDM3Y, the effective nucleon-nucleon interaction v{r) is assumed to be density and energy dependent 
and therefore becomes functions of density and energy and is given by 

v{r,p,E)=t'^^^{r,E)g{p,E) (4) 

where t^^^ is the same M3Y interaction given by Eqn. 3 but supplemented by a zero range pseudo-potential [18] 

^M3Y ^ 79996 _ 2134!^ + Jm{E)5{r) (5) 
4r 2. or 

where the zero-range pseudo-potential representing the single-nucleon exchange term is given by 

jQf^{E) = -276(1 - 0.005£;/A)MeV.fm^ (6) 
and the density dependent part has been taken to be [19] 

g{p,E) = C{l-(}{E)p'h (7) 

which takes care of the higher order exchange effects and the Pauli blocking effects. E/A is energy per nucleon. The 
constants of this interaction C and [3 when used in single folding model description, can be determined by nuclear 
matter calculations [20] as 2.07 and 1.624 fm^ respectively. The density dependence parameter (i has the dimension of 
cross-section. The term (1 — I3p'^/^) reduces the strength of the interaction and changes sign at high densities making 
it repulsive. This is a direct consequence of the Pauli blocking effect. Thus (1 — /3p^/^) represents the probability of 
non-interaction arising due to collision probability fif?!"^ of a nucleon in nuclear medium of density p. The parameter 
/3 can be identified as the 'in medium' effective nucleon-nucleon interaction cross section oq. This value of /3 along 
with nucleonic density of infinite nuclear matter po can also provide the nuclear mean free path A = (poco)~^- 

III. CALCULATION AND ANALYSIS 

The nuclear ground state densities have been calculated in the framework of spherical Hartree Fock plus BCS 
calculations in co-ordinate space using two different parameter sets of the Skyrme interactions given by [21] 

Wrai(ri, ra) = to(l + xqP ^)b{v) + ti(l + XrV^)\b(v)\? + k'^(5(r)] + 1^(1 + X2P.)k'(5(r)k + ^3(1 + a;3Pa)p"(R)'5(r) 

(8) 

There are negligible differences in the ground state densities of the nuclei with the two different parameter sets (for 
SIII [22] and SkM* [21]) provided in Table 2. We show SkM* parameterization for the ^^O ground state density (Fig. 
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1) and used it for calculations of the single folding potential and form factor. We carried out the same procedure for 
""^^Ne and for SkM* the binding energies per nuclcon arc in good agreement (within 0.1 MeV) with the experimental 
values for both the nuclei. 

The potentials used while calculating phenomenological best fits [6,23] have the following form, 

l^pheno(r) = -Vo fo{r) - i Mr) + 4 i a,Ws{d/dr)fs{r) + 2{n/m^cf Vs.ol/r{d/dr)Uo{r) (L.S) + Koui (9) 

where, fxir) = [1 + exp C"'^" )]~-^, = VxA^I"^ and x = o,v,s,s.o. The subscripts o,v,s,s.o denote real, 
volume imaginary, surface imaginary and spin-orbit respectively and Vo, Wv {Wg) and Vg.o are the strengths of the 
real, volume (surface) imaginary and spin-orbit potentials respectively. K;oui is the Coulomb potential of a uniformly 
charged sphere of radius 1.20 A^/^. 

In semi-microscopic analysis both the volume real {V) and volume imaginary (W) parts of the potentials (generated 
microscopically by folding model) are assumed to have the same shape, as in Ref. [4], i.e. Knicro(?') — V + iW = (TVr, 
-|- iNi)U(ri) where, Nji and Ni are the renormalization factors for real and imaginary parts respectively [9]. In Fig. 
2 the p + renormalized real folded potentials (employing the SBM and DDM3Y interactions) are shown at E 
= 24:.5A MeV, along with the best fit phenomenological real part. Thus the potentials for elastic scattering analysis 
include real and volume imaginary terms (folded potentials) and also surface imaginary and spin-orbit terms (best fit 
phenomenological potentials). 

r n 2 

For each angular distribution, best fits are obtained by minimizing x^/^j where = Sfe=i ""'^'^aI '(eT)^^''^ ' 
where ath and aex are the theoretical and experimental cross sections respectively, at angle 0k, Aaex is the experimental 
error and N is the; nmnber of data points. The potentials for elastic scattering analysis are subsequently used in the 
DWBA calculations of inelastic scattering with transferred angular momentum /. The calculations are performed 
using the code DWUCK4 [24]. The derivative of the potentials (S^) are used as the form factors. The microscopic 
real and imaginary form factors have the same shape with strengths and Nf^ respectively, where A^p^ = 

^R,ifYms^ where the radius parameter is the rms radius of the folded potential. In addition, form factors derived 
from phenomenological surface imaginary and spin-orbit potentials are included. The deformation parameters S are 
determined by fitting the inelastic scattering angular distribution. The renormalizations required for the potentials 
are reminiscent of those for deuteron and ^Li scattering and it may be surmised that it is for the same reasons; 
weak binding and ease of breakup and other reaction channels. Table 1 gives the parameters of the interactions used 
here. Both interactions provide incompressibility of ^ 300 MeV for spin and isospin symmetric cold infinite nuclear 
matter. Moreover, in case of DDM3Y, (3 value of 1.624 fm^ obtained from nuclear matter calculations, along with 
nucleonic density of 0.16 fm~^ provides a mean free path of about 4 fm which is in excellent agreement with other 
theoretical estimates [25]. Table 3 gives the phenomenological best fit optical model parameters while Table 4 gives 
the renormalization factors, S values, x^/N for the folding model analysis. 

The relationship between the reduced electric quadrupole transition rate B{E2) for the ground state to the 2+ state 
in units of fm^ and the quadrupole deformation parameter 5 is given by [26] 

^(1 + 0.16^ + 0.20(5^ + ....) = A'KB{E2f''^/{iZE?) (10) 

where R = 1.2^^/^ fm and Z is the atomic number. The quadrupole deformations listed in reference [26] were obtained 
by using Eqn. 10 but keeping only the terms up to first order in S. We have recalculated these values by keeping 
terms up to third order in 6. The recalculated quadrupole deformations thus obtained from the experimental B{E2) 
values listed in reference [26] are 0.33078 and 0.59284 for and ^^Ne respectively. As can be seen from Table 4, 
the quadrupole deformation obtained from the present analysis for ^^O is in excellent agreement while that for ^^Ne 
is significantly underestimated due to lack of experimental data at forward angles. The inelastic scattering is more 
sensitive at forward angles due to its relative purity compared to data corresponding to backward angular range where 
other non-elastic processes also contribute. The quality of elastic as well as the inelastic fits deteriorate at backward 
angles and similar deterioration of fits are also seen for proton scattering from other nuclei [10,27]. The reason for this 
is probably that the full cross section is ascribed to potential scattering while quasi-compound- nucleus formation feeds 
back into the elastic channel and whose energy-dependence is controlled by barrier-top resonances. The backward 
angular range classically corresponds to smaller impact parameters. This fact suggests higher compound nuclear 
formation probabilities at backward angles while those at forward angles are expected to be negligibly small. Since 
relative contributions of compound elastic and direct elastic arc not disentangled, x^ was calculated only upto about 
90° in center of mass during fitting the ^^0 data. 

It may be noted that since we could not acquire the experimental data for ^*0, they have been read quite accurately 
from the plots in the original paper [6] and a 5% uniform error has been assumed. The technique used for extracting 
the data from plots along with relevant co-ordinate transformations are described in the appendix. 
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IV. CONCLUSION 



In the present study we find that the parametrized SBM effective interaction, and the reahstic DDM3Y effective 

interaction obtained from sophisticated G-matrix calculations provide equivalent descriptions for the clastic and 
inelastic scattering of protons from the mirror nuclei ^*Ne and ^^O (Fig. 3, 4). The values of the deformation 
parameters have been extracted from the calculations for these nuclei. Even though the analysis reported here is 
quite detailed, measurements at higher energies may still be useful in distinguishing various effective interactions. 
The lack of sensitivity is due to the fact that at such energies, nuclear densities primarily probed are near the surface, 
which makes the density dependent effects less realizable. The form of density dependence used here is more physical 
as compared to other forms. The parameter (3 can be interpreted as 'in medium' nucleon-nucleon interaction cross 
section while (1 — (ip^/^) as the non-interaction probability arising due to higher order exchange and Pauli blocking 
effects. 

In summary, a consistent folding model analysis of proton scattering on A = 18 nuclei is carried out using two 
different effective NN interactions. The conventional way of generating the form factors is followed, that is, by taking 
the derivatives of the potentials (microscopic real and imaginary as well as phcnomcnological surface imaginary and 
spin-orbit potentials). Deformation parameters {6) are extracted from the analyses. The results obtained for the 
deformation parameter are in good agreement with the available results. The density dependence parameter obtained 
from nuclear matter calculations, which has been used in the single folding model description for the analysis of elastic 
and inelastic scattering of protons, also provides excellent estimate for nuclear mean free path. We know that the 
DDM3Y effective interaction, has profound theoretical basis. It provides unified description of cluster radioactivity, 
scatterings of a and heavy ion when used in a double folding model, and nuclear matter when used in a single folding 
model. We find that it also provides reasonable description for elastic and inelastic scattering of protons. 

The authors gratefully acknowledge L. A. Riley for sending the experimental data in a tabular form. 



Table 1: 

Parameters of the SBM and DDM3Y interactions 





SBM 


n 


Ci(MeV) 


C„(McV) 




a(fm) 


6(MeV/c) 




d(fm) 




2/3 


215.7 


669.3 




0.554 




668.7 






0.813 


DDM3Y 


n 




C 


a(MeV-' 


) 


/3(fm^) 














2/3 


2.07 


0.005 




1.624 














Table 2: 
























Parameter sets of the Skyrme interactions 




















Interaction 


to 


ti 








Xo 


Xl 


X2 






Q 


SIII 


-1128.75 


395.0 


-95.0 


14000.0 


120.0 


0.45 


0.0 


0.0 


1.0 




1 


SkM* 


-2645.0 


410.0 


-135.0 


15595.0 


130.0 


0.09 


0.0 


0.0 


0.0 
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Table 3: 

Phenomenological potential parameters used in p -|- ^^Ne and ^^O scattering 



Nucleus 


E/A 


Vo 


To 


fflo 




rv 




Ws 




as 


Vs.o 




ffls.o 


Xci/N 


J/A 


Ref 




(McV) 


(McV) 


(fm) 


(fm) 


(McV) 


(fm) 


(fm) 


(McV) 


(fm) 


(fm) 


(McV) 


(fm) 


(fm) 




(MeV fm^) 






30.0 


40.00 


1.100 


0.730 








7.00 


1.380 


0.600 


7.80 


1.090 


0.740 


5.265 


-364.7 


[6] 


18Q 


24.5 


48.57 


1.163 


0.780 


2.266 


1.169 


0.690 


5.77 


1.169 


0.690 


5.90 


0.882 


0.630 


6.896 


-527.0 


[23] 
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Table 4: 

Renormalizations of SBM and DDM3Y folded potentials and form factors for p + ^*Ne and scattering at incident 
energy [E/A) and excited state energy (E*) in MeV, angular momentum transfer (I), deformation parameter (6), 
volume integral {J/A) of the real folded potential in MeV fm^ and x^/N values from best-fits to the elastic and 
inelastic scattering data 





Nucleus 


E/A 


E* 


^R 


iVi 


' rms 






/ 


5 


Xci/N 


X?nel/N 


J/A 


Interaction 




30.0 


1.890 


0.67 


0.000 


3.462 


2.320 


0.000 


2 


0.436 


5.992 


0.973 


-336.1 


SBM 




30.0 


1.890 


0.62 


0.000 


3.538 


2.194 


0.000 


2 


0.400 


5.798 


0.859 


-358.0 


DDM3Y 


i8Q* 


24.5 


1.982 


0.87 


0.006 


3.547 


3.086 


0.020 


2 


0.400 


7.131 


4.471 


-502.7 


SBM 


18q* 


24.5 


1.982 


0.80 


0.006 


3.508 


2.806 


0.019 


2 


0.367 


6.704 


7.696 


-469.2 


DDM3Y 
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V. APPENDIX 



This appendix is aimed at providing a simple co-ordinate transformation formula for general purpose use in ex- 
tracting data from plots using computer graphics. 

In many instances the values of the experimental data are not available in literatures while the plots are shown. 
Particularly, for old plots, tabulated data are often very difficult to acquire. In such cases, data values can be read 
in quite accurately by using computer programs. In some cases both linear and logarithmic scales are involved. A 
general formula is thus required which will convert the co-ordinates of the data points as shown by the software to 
their actual values. It should also be taken into account that the co-ordinate system of the printed graph may be 
rotated, albeit small, with respect to the software co-ordinate system. Thus a handy formula would be extremely 
useful for such purposes. 

Wc assume {X, Y) to be the co-ordinates of the printed graph while {x, y) arc those of the computer program such 
as GSview. In Fig. 5, {xQ,yo) is the origin of the printed graph while (xi,yi) and {x2,y2) are two points on the X 
and Y axes, respectively. 

Since the ratio of two linear lengths measured in one co-ordinate system should be equal to that in the other system, 

X-Xq _ {x-xo)sece+[{y-yo)-{x-xo)tan0]sin0 Vi-yo . i-^-^-i • i +• 

y _y — ^ f — _ \ a — where, tanO = _ ■ A little manipulation gives, 

Y V ^iv V \ \ ^^~ xo)cosO + {y- yo)sin0 ^ 

X = Xo + {X,- Xo) [ _ ^,),ec^ \ (11) 

We consider X scale to be linear and Y scale to be logarithmic. Therefore, 

logY-loqYo ^ \{y-yo)-{x-xo)tan0]cose ^^^^^ 
logY2-logYo {y2-yo)sec0 ^ ' 

-(x-j;o)si.ie+(y-yo)coa9 

Y = Yo(g-) (12) 

If both axes are linear, then Eqn. 11 holds while 

Y V j.(v ^>,\ -ix-xo)sin0+{y -yo)cos0 i^ 

Y — Yo + [^2 - Yo) -. r 

L (y2 - yo)sec0 J 

If both axes are logarithmic then Eqn. 12 holds while, 

X = Xo(|i) (14) 

The expressions provided above, though simple, will be immensely useful for experimentalists as well as theoreticians 
for extracting data where their explicit values are not available. 
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FIG. 1. Skyrme Hartree-Fock densities of O used in this work. 
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p + E = 24.5 MeV 
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SBM 
DDM3Y 



FIG. 2. The phenomenological and renormalized real folded potentials (employing SBM and DDM3Y interactions) for p + 
O at 24.5A MeV. 
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FIG. 3. The experimental angular distributions and folding model calculations (SBM) of (a) p + '^^Ne at 30A MeV for 

elastic and inelastic [E* = 1.890 McV (2+)] scattering [5], (b) p + ^**0 at 24.5A McV for clastic and inelastic [E" = 1.982 
MeV (2+)] scattering [6]. The corresponding A^r, Ni, , values and phenomenological surface imaginary and spin-orbit 



parameters are given in Tables 3, 4. 
sections respectively. 



The continuous and dashed lines correspond to calculations for elastic and inelastic cross 
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FIG. 4. The same as in Fig. 3 but employing DDM3Y interaction. 
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FIG. 5. {X, Y) are the printed graph co-ordinate system, rotated by an angle 6 with respect to {x, y), the computer program 
co-ordinate system. 
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